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Abstract — We propose a new code design that aims to distribute 
an LDPC code over a relay channel. It is based on a split- 
and-extend approach, which allows the relay to split the set 
of bits connected to some parity-check of the LDPC code into 
two or several subsets. Subsequently, the sums of bits within 
each subset are used in a repeat-accumulate manner in order to 
generate extra bits sent from the relay toward the destination. We 
show that the proposed design yields LDPC codes with enhanced 
correction capacity and can be advantageously applied to existing 
codes, which allows for addressing cooperation issues for evolving 
standards. Finally, we derive density evolution equations for the 
proposed design, and we show that Split-Extended LDPC codes 
can approach very closely the capacity of the Gaussian relay 
channel. 

I. Introduction 

By exploiting the broadcast nature and the inherent spa- 
tial diversity of wireless communications, Sendonaris et al. 
introduced the concept of cooperative diversity IH, Q over 
wireless relay channels and their multi-terminal extensions. 
A relay channel is a three terminal network consisting of 
a source, a relay, and a destination. The source broadcasts 
a message to both relay and destination, while the relay 
forwards the message to the destination. Subsequently, many 
authors proposed cooperation protocols for the relay channel, 
which can be classified into two major categories, namely the 
amplify-and-forward (AF) and the decode-and-forward (DF) 
O. In AF protocols, the relay simply amplifies the received 
signal and forwards it to the destination. The DF protocol 
allows the relay to decode the received signal, re-encode it, and 
forward it to the destination. The forwarded message can either 
be identical to, or part of the initial transmission (repetition 
coding), or it can be obtained by using a dedicated coding 
scheme at the relay (distributed coding). In the first case the 
destination combines received signals from both source and 
relay, which results in an improved signal-to-noise ratio (SNR) 
on the received transmission. Besides, the same code is used 
for encoding at the source and decoding at the destination. 
In the second case, the destination gains knowledge of extra 
information, but it needs a dedicated decoding scheme, able 
to jointly decode received signals from both source and relay. 

One of the most known examples of distributed coding is 
the one of a distributed turbo-code [4]: the source broadcasts 
a recursive convolutional code (RCC) to both relay and des- 
tination. After decoding, the relay interleaves and re-encodes 
the message using the same RCC, prior to forwarding it to 
the destination. Because the destination receives both codes 
in parallel, it can jointly decode received signals from source 
and relay by using a parallel-concatenated turbo-code. 

This work was partially carried out in the scope of Celtic CP5-026 
(WINNER+) project. 



Low-Density Parity-Check (LDPC) codes play a prominent 
role in the family of error-correcting codes. They feature low 
complexity decoding and can be optimized for a broad class 
of channels, with performance approaching the theoretical 
Shannon limit [5J. Although they lend themselves less easily 
to distributed schemes, several approaches have been already 
proposed in the Hterature fTl, fSl, [51, [TOl, {\T\, {VT\, 
[13|. Some of these approaches are somehow based either on 
serial or parallel code concatenatioifl or on punctured (rate- 
compatible) LDPC codes. From the code design point of view, 
the serial or parallel concatenation of LDPC codes has intrinsic 
limitations, mainly because parity-check matrices used for 
decoding at the relay and at the destination are included one in 
the other, resulting in inappropriate matrix topologies (density 
on non-zero entries, column and row weight distributions, 
cycles, etc.). Punctured codes also present some weaknesses 
in the context of coded cooperation. In such a cooperation 
scheme, punctured bits, which are not broadcasted by the 
source, are transmitted from the relay toward the destination. 
Hence, the punctured code has to be robust, such as to allow 
successful decoding at the relay: in practice punctured bits 
are those bits which are the easiest to retrieve by the iterative 
decoding process. On the other hand, bits transmitted from the 
relay to the destination encounter better channel conditions: 
advantageously, these bits should be those which are the 
most difficult to retrieve by the iterative decoding process. 
The contradiction between these two requirements leads to 
an imbalance in the design of the puncturing pattern. 

In this paper we propose a new code design method that 
aims to create incremental redundancy for LDPC codes, while 
avoiding both code concatenation and code puncturing. It is 
based on a split-and- extend approach, which can be seen as the 
"coding analogous of the divide-and-conquer concept". After 
decoding the received signal, the relay computes extra parity 
bits by splitting parity-checks of the initial code. Hence, each 
extra parity bit is the sum of some subset of bits participating 
in the same parity-check of the initial code. Then the relay 
transmits these new parity bits towards the destination. The 
whole process amounts to create a new matrix, whose rows 
correspond to parity-checks involving both old and new parity 
bits. These parity-checks are therefore distributed over the 
relay channel, in the sense that part of checked bits are 
received on the source-to-destination link, and another part 
are received on the relay-to-destination link. Consequently, this 
new matrix can be used at the destination to jointly decode 
the received signals from both the source and relay. 

'Meaning that the graph of the LDPC code broadcasted from the source is 
a subgraph of the destination decoding graph. 



The paper is organized as follows. The proposed code 
design is introduced in Section In Section |III] we show 
that the proposed design can be advantageously applied to 
to existing codes, which allows for backward compatibility 
while addressing cooperation issues for evolving standards. 
In Section |IV] we propose a "coding-perspective" analysis 
of cooperative systems. We introduce threshold and capacity 
functions, and we derive density evolution equations for Split- 
Extended (SE) LDPC codes over the Gaussian relay channel. 

Numerical results are shown in Section |V] and Section |VI] 
concludes this paper 

II. Split-Extended LDPC codes 
A. Basic idea 

The basic idea of the split-extend design can be resumed as 
follows. Let Hi be the parity-check matrix of the LDPC code 
broadcasted by the source to both relay and destination. Hence, 
broadcasted bits satisfy parity-check equations corresponding 
to the rows of Hi . After decoding the received signal, the relay 
computes extra parity bits by splitting these parity-checks, 
as illustrated at Figure [T] The parity-check in the middle 
corresponds to a row of Hi. In the left example, a new parity 
bit ei is created by spitting the original parity-check into two 
sub-checks. Precisely, this means that the set of bits connected 
to the check-node is partitioned into two subsets, and the parity 
bit ei is generated as the sum of the bits of either one of the 
two subsets. In the right example, two new parity bits ei and 
62 are created by spitting the original parity-check into three 
parity-checks. Precisely, the set of bits connected to the check- 
node is partitioned into three subsets, and ei is generated as 
the sum of the bits in the first subset. Subsequently, e2 can be 
generated either as the sum of ei and the bits in the second 
subset, or as the sum of the bits in the third subset. The total 
number of extra parity bits depends on the number of rows of 
Hi and the number of extra bits generated for each row of Hi 
(which may vary from one row to another). The sequence of 
all the extra parity bits, denoted by E ^ (ei, 62, . . . ), is then 
transmitted from the relay to the destination. The matrix H 
obtained by the split-extension of Hi (i.e. the incidence matrix 
of the split-extended graph) verifies H -{X, EY =0, where X 
denotes the codeword broadcasted by the source. Therefore, 
H can be used at the destination in order to jointly decode 
the received signals from both source and relay. 

A more general example of split-extension is illustrated at 
Figure |2] The original check-node is split into several sub- 
checks, and extended bits are generated in a repeat-accumulate 
manner Such a split-extension will be referred to latter in the 
paper (Section llVb as repeat-accumulate split-extension. 




Fig. 1. Split-Extension examples 



B. The general case 

The following notation wiU be used throughout this section: 

• For any positive integer A^, [1 : TV] = {1, . . . , denotes 
the set of integers between 1 and N, inclusive. 

• For any subset 5* C [1 : iV], [1 : iV] \ 5 denotes the set of 
integers between 1 and N that are not in S. 

• For any length-iV vector V — . . . , ujv) and any subset 
S = {ii,...,ik} C [1 : iV], V\s denotes the length-A: 
vector defined by the coordinates of V which are in S, that 
is V\s = {v^,,...,Vi^^) 

• For any binary matrix H of size M x N, TZ{H) C {0, 1}^ 
denotes the set of the M row-vectors of H. A partition of 
TZ{H) is a set of nonempty subsets Vi, . . . , Pmi C TZ{H) 
such that every row of H is in exactly one of these subsets. 

• For any binary vectors Vi, V2 G {0, 1}^, Vi © V2 denotes 
their componentwise sum modulo 2. 

[Definition] Let Hi and H be two parity-check matrices of 
size Ml X A^i and AI x N, respectively, with A^i < A^. We say 
that the matrix H is obtained by split-extending the matrix Hi 
if there exist a partition of the M rows of H in Mi disjoint 
subsets Pi, ... , Vmi ^ 7^(i^), and a subset 5 C [1 : A^] with 
cardinality A^i, such that for any to G [1 : Mi]: 

i|s = i?mand L|[i:7vi\s = 0, 

where Ri, . . . , Rmi denote the Mi rows of Hi. In this case, 
the s,tt E — [\ : N] \ S is called set of extended bit-nodes. 
The split-extension is called non-singular if the columns of H 
corresponding to E are linearly independent. 

Now, consider two parity-check matrices Hi and H, such 
that iJ is a non-singular split-extension of Hi. Then: 

• If X is a codeworc0 of H, then X\s is a codeword of Hi 

• For any codeword Xi of Hi, there exists a unique codeword 
X of H, such that X\s ^ Xi 

Matrices {Hi, H) can be used within a cooperative trans- 
mission system as follows: 

• The source encodes the packet of information bits, gener- 
ating a codeword Xi of Hi . It broadcasts Xi to both relay 
and destination. 

• The relay decodes the received signal, correcting the trans- 
mission errors on Xi. It generates a codeword X of H, 
such that X\s = Xi, and sends the set of extended bits 
X\e towards the destination. 

• Thus, the destination receives noisy versions of X\s = Xi 
and X\e (from both the source and the relay), which can 
be decoded using the matrix H. 

-By abusing language, we say that X is a codeword of H, if HX^ = 




Fig. 2. Repeat-Accumulate Split-Extension 



III. Split-extend design for backward 

COMPATIBILITY 

This section is independent of the following sections, 
though, it highlights an interesting property of the proposed 
design: it can be advantageously applied to existing codes, 
which allows for addressing cooperation issues for evolving 
standards, while maintaining backward compatibility with a 
reduced impact on user equipment. To illustrate this, the LDPC 
codes from the IEEE.802.16e (WiMAX) standard CH with 
coding rates 1/2 and 2/3 have been split-extended, such that 
the number of generated extended bits be equal to the number 
of information bits. Thus, for coding rate 1/2, each row of 
the parity-check matrix has been split into two rows; while for 
rate 2/3, each row of the parity-check matrix has been split 
into three rows (see Figure [U. Splitting has been performed 
by a dedicated algorithm that search for short cycles in the 
parity check matrix, then splits rows such that to break as 
many short cycles as possible. Base matrix of the Quasi-Cyclic 
(QC) LDPC WiMAX code with rate 1/2 and the corresponding 
split-extended matrix are shown at Figure |3] (—I's entries of 
the base matrix are represented by a dash sign). 

Clearly, split-extended matrices can be used to address 
cooperation issues for uplink transmissions, in a completely 
transparent way for the user: the user encodes the transmitted 
signal by using the original parity check matrix; the relay 
decodes the signal, then computes and sends the sequence of 
extended parity bits to the base station, which will use the 
split-extended matrix in order to decode the received signals 
from both user and relay. Hence, split-extended matrix is only 
needed at the relay and the base station. For downlink trans- 
missions, if the user terminal is equipped with split-extended 
matrices, the situation is symmetric. Otherwise, the relay can 
only repeat the sequence of information bits, which, however, 
provides the user with an energy gain on the information 
sequence. 

Simulation results over the AWGN relay channel, with 
QPSK modulation, are shown at Figures |4] and |5] The source 
broadcasts either the WiMAX code with rate 1/2 (Fig. [Hi or 
the WiMAX code with rate 2/3 (Fig. EJ, and the SNR on the 
source-to-relay link is fixed to 2.5 and 4.5 dB, respectively. 
Both figures compare the performance of two cooperation 
scenarios: the relay generates and transmits extended bits 
in the first scenario, while in the second, it forwards the 
(error-corrected) sequence of information bits. Plotted curves 
represent SNRs required on source-to-destination and relay-to- 
destination links, such that to obtain a target frame error rate 
either of 10^^ (dotted curves) or of 10^^ (solid curves). The 
Self-Corrected Min-Sum algorithm |15| is used for decoding 
at both relay and destination. The gap between the dotted and 
solid curves is determined by the slope of the frame error 
rate curves in the waterfall region. However, for the second 
scenario (repetition of the information sequence) the increased 
gap between the two curves for small SNRsd values is also 
justified by a frame error rate error floor above lO"**. Note 
also that only SNR pairs with SNRrd > SNRsd are likely to 
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Fig. 3. Base matrix of the QC-LDPC WiMAX code with rate 1/2 (top), and 
corresponding split-extended base matrix (bottom) 
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be encountered in practice. The SNR gain between the two 
scenarios can be measured either as the horizontal distance 
(for the source-to-destination link) or the vertical distance (for 
the relay-to-destination link) between corresponding curves. It 
can be observed that split-extended codes achieve a significant 
SNR gain, in order of several dBs, over the repetition scenario. 



IV. Asymptotic analysis of SE-LDPC codes 

We denote by £{X,p) the ensemble of LDPC codes with 
edge-perspective degree distribution polynomials A and p lfT6l . 
It is well known that when the code length tends to infinity, 
(almost) all the codes of the family behave alike, and they 
exhibit a threshold phenomenon, separating the region where 
reliable transmission is possible from that where it is not |fT6l . 

Consider some channel model depending on a parameter 
(T, such that the channel conditions worsen when a increases 
(for instance, the noise variance for the AWGN channel, or 
the error probability for the BSC channel). The threshold 
of the ensemble £{X, p) is defined as the supremum value of 
a (worst channel condition) that allows transmission with an 
arbitrary small error probability, assuming that the transmitted 
data is encoded with an arbitrary-length code of £{X,p). 

The above threshold phenomenon can be extended to the 
relay channel, but the following must be taken into account: 

• the channel is modeled by three parameters crsR. ctsd- and 
(Trd, corresponding, with obvious notation, to the three links 
between source, relay, and destination. 

• the ensemble of SE-LDPC codes depends not only on A and 
p, but also on the splitting distribution. 

[Channel assumptions] Since we are strictly interested on code 
analysis, the following assumptions will be made: 

• when the relay fails to decode the received signal from 
the source, it does not transmit any information to the 
destination, 

• the relay channel is degraded, in the sense that the above 
parameters must satisfy ctsr < ctsd and ctrd < ctsd- 

[Distributed code] A linear distributed code of dimension K 
is a vector subspace C C F^^ x F^^, such that C and its 
projection on F^^ are both of dimension K. The distributed 
rate of C is by definition (ri, — ■^)- Hence, ri < 1, 
but r2 can be greater than 1. The overall coding rate is defined 
as r = .r ^ — ^'^P . The idea behind is that the first A^i 

-'V1+-1V2 ri+r2 ^ 

bits of a codeword c e C are broadcasted from the source to 
both relay and destination and, in case that the relay manages 
to decode the received signaj^ it transmits the last iV2 bits 
toward the destination. 

[SE-LDPC ensembles] Let Hi be the parity-check matrix of 
the LDPC code broadcasted by the source, and let 9 > 2. For 
each parity-check of Hi, assume that: 

• the set of bits connected to the parity-check is partitioned 
into 9 subsets of (almost) regular size; that is, each subset 
contains either [|J or [|] bits participating in the parity- 
check, where d denotes the parity-check degree, 

• extended bits are generated in a repeat-accumulate manner 
(Figure |2]). 

The resulting distributed SE-LDPC code will be referred to as 
having splitting degree 9. We denote by £g{X, p) the ensemble 
of repeat-accumulate SE-LDPC with splitting degree 9, ob- 
tained by split-extending LDPC codes with edge-perspective 
degree distribution polynomials A and p. 

^Thus, this definition is dependent on the above channel assumptions. 



A. SE-LDPC ensemble thresholds 

Given an ensemble of SE-LDPC codes, our intention is 
to separate the region of parameters (ctsr, ctsd, ctrd) where 
reliable transmission is possible from that where it is not. 

Let al{X,p) be the threshold of the £{X,p) ensemble. If 
fSR > the error decoding probability at the relay is 

lower bounded by a positive constant. Consequently, reliable 
cooperation cannot be achieved, as the relay does not transmit 
any information to the destination when it fails to decode the 
received signal, and the destination cannot reliably decode the 
signal received from the source, since (Tsd > ^sr > crJ(A, p). 

From now on we consider that ctsr < (tI{X,p); hence, we 
may assume error freJ^ transmission between source and relay. 
We will also use the following notation: 
9 (J — (jRD (we drop subscript rd), which will be referred to 

as noise parameter, 

• 5 = > 1, which will be referred to as channel 

ctrd 
discrepancy. 

The threshold function (Jg x p ■ [1 1 00 [— > M associates with 
each discrepancy value 6 the noise threshold ag ^ ^{6), defined 
as the supremum value of a that allows transmission with an 
arbitrary small error probability, assuming that the transmitted 
data is encoded with an arbitrary-length distributed code from 
£e{X,p). This definition makes sense only under the implicit 
assumption of a concentration theorem as in [161, which can 
indeed be derived by using the same arguments as in loc. cit. 

The threshold function can be efficiently computed by 
tracking the density evolution of messages exchanged within 
the iterative decoding, as explained in the next section. 

B. Density evolution 

Throughout this section, we assume binary-input AWGN 
relay channel. We combine a multi-edge approach ifTTl and 
the Gaussian approximation method proposed in [18), in order 
to derive density evolution equations for the SE-LDPC code 
ensemble £e{X,p). We separate the set of bit-nodes of the 
expanded graph into two subsets: 

• type-1 bit-nodes, which correspond to bits received by the 
destination from the source, 

• type-2 bit-nodes, which correspond to extended bits received 
by the destination from the relay. 

We distinguish between type-1 and type-2 edges, according 
to whether they are incident to type-1 or type-2 bit-nodes. 
Moreover, check-node degrees are also defined type-wise. 
Hence, we say that a check-node is of degree {di,d2) if it 
is connected to di type-1 bit-nodes and d2 type-2 bit-nodes. 
From our definition of SE-LDPC ensembles, it follows that 
the type-2 degree c?2 is equal either to 1 or 2 (see Fig. |2]). 
Finally, for each type i = 1, 2, we define: 

• a|j*' is the fraction of type-t edges connected to bit-nodes 
of degree d, 

• is the fraction of type-t edges connected to check- 
nodes of degree (c?i,d2). 

*With arbitraiily small eiTor probability, as the code length tends to infinity. 
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^ is the average check-node degree of the 
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C 0, if 2 < j < 6* - 1 
kj,x = I I, if j = 1 or j = + 1 
2, if - 6» + 1< 7 < 



kj,2 = 0- \j\ - kj,i 



Now, under the Belief-Propagation decoding, let Ti^jt] de- 
note the mean of outgoing messages from type-i bit-nodes 

It] f \ 

at iteration I. Let also ry = 1 — E [ tanh — ] , where E 

denotes as usual the expected value operator. Define: 




d«, ((/){0) = 1) 



11,22 
[2] 



Then, under the assumption that the messages exchanged 
during the iterative Belief-Propagation decoding are indepen- 
dent and symmetric Gaussian distributed, rf^ and rf^ can be 
recursively computed by: 



with initial values 



.[11 J-2\ 



The 



proof will be omitted, since it follows from the same argu- 
ments as in ifTSl . The above recursion holds as long as i is 
less than half the girth of the graph, which goes to infinity 
with the code-length, and the successful decoding condition 
for an "infinite" code from £e{X,p) can be expressed as 
('"i '> '^f ) ~^ 0- Therefore, the threshold function defined in 
the above section, can be computed by: 

o-e,A,p('5) = sup{CT I lim = lim = 0} 

C. Theoretical limit 

In order to evaluate the performance of an ensemble of 
codes, we would like to compare its threshold function with 
the capacity function, inferred from the channel capacity. 

Capacities of various relaying strategies have been com- 
puted in [|19l , II20I . II2TI . and depend on the capacities of the 



three links. Since we assumed that source-to-relay transmis- 
sion IS error frell we only consider the two other links. Let 
7rd(<7) denote the information rate capacity of the relay-to- 
destination channel with parameter a, and let 7sd(o') be de- 
fined in a similar manner The information rates are considered 
by transmitted bit, thus both 7rd(c), 7sd(o') £ [0, 1]. We also 
assume that the noise parameter a G [0, +oo[ and 7rd, 7sd are 
continuous decreasing functions, such that 7rd(0) = 7sd(0) = 
1 and lim 7rd(o') = lim 7sd(c) = 0. 

Now, let (ri,r2) € [0, 1] x [0, +oo[ be a target distributed 
coding rate. The capacity function 7ri,r2 ■ [I7 +00 [0, +oo[ 
is defined by 7ri,r2(<^) — '^^ where a is the unique solution of 
the equation: 

7sd(<5o-) 7rd(o-) 



= 1 



r2 



Note that for cr = 0, we have 



> 



7Sd(0) I 7RD(n) ^ J_ I J. 

r2 ri r2 

i > 1, and lim 'Jsd!M + ImiZl = 0, thus such a solution 

■'•1 - '■2 

always exists and it is unique, due to the above assumptions. 

The meaning of the capacity function is the following. 
Assume that we want to transmit information with distributed 
rate (ri,r2) over some relay channel. The rate ri is chosen 
according to the quality of the channel between source and 
relay, such that to ensure error free transmission between 
them. The rate r2 is generally chosen according to the delay 
constraints of the cooperation system. The question is to know 
is there exists a distributed code with distributed rate (ri , r2 ) 
allowing error fre^ transmission. The answer is as follows. 
If the discrepancy and noise parameters {S, a) of the relay 
channel verify a < jri.r2 {^) then such a distributed code exits. 
However, note that a code allowing error free transmission for 
some pair {S, a), might not be suitable for some other pair of 
parameters satisfying the above condition. If > 7ri,r2(^)5 
then reliable transmission with distributed rate (ri , r2 ) is not 
possible. The proof will be given in an extended version of 
this paper 

V. Numerical results 

We assume BI-AWGN relay channel throughout this sec- 
tion. The following degree distribution pair, with designed 
coding rate 1/2, was designed by exact density evolution, and 
its thresholcfl over the BI-AWGN channel is a* ^ 0.9649 [5 |. 

\{x) = 0.2199j: + 0.2333x^ + 0.0206j:^ + 0.0854a;^ + 0.0654x'^ 

+0.0477i'' + 0.01912;** 4- 0.08062:^** + 0.22802:1^ 
p{x) = 0.6485a:'^ + 0.3475a;« + 0.0040x^ 

First, we consider the SE-LDPC code ensemble £2{X,p). 
The designed distributed rate is (ri,r2) = (1/2, 1), meaning 
that the source broadcasts a code with rate 1/2, and the relay 
generates a number of extended-bits equal to the number 
of information bits. If the standard deviation of the white 
Gaussian noise on the source-to-relay link is less than the 

^In practice, the channel need not be error free; the assumption is that the 
(Tg]^ noise is below the threshold of the code broadcasted by the source. 
^Arbitrary small error probability when the code length goes to infinity. 
'Note that the threshold calculated by Gaussian approxiination is 0.9459. 



above threshold a* or, equivalently, the signal to noise ratio 
is greater than SNR* = —2.70 dlio we can assume error free 
transmission between source and relay. From our definition 
of the discrepancy, the signal to noise ratios on the source- 
to-destination and relay-to-destination links, are related by 
SNRrd = SNRsD + A, where A = lOlogiol'?^) is the 
discrepancy value in dB. The ensemble threshold and the 
capacity functions are plotted at Figure |6] We can observe that 
the gap between the two curves is relatively small (between 0.7 
and 0.3 dB) for discrepancy values A e [0, 8.5], and it begins 
to increase starting from this point. Hence, if a discrepancy 
value A > 8.5 dB is not likely to be encountered in practice, 
the above SE-LDPC code can be used to achieve reliable 
communication for channel parameters (5, a) very close to 
the capacity. 

[Remark] HARQ systems with incremental redundancy rep- 
resent another possible application of the proposed SE-LDPC 
codes. In this case extended-bits are transmitted by the source 
as incremental redundancy, whenever the destination fails to 
decode the originally received signal. In such a case, the 
discrepancy is expected to take on relatively small values. 

Figure |6] shows also the threshold function for the SE- 
LDPC code ensemble Ez(\tP\ whose designed distributed 
rate is (ri,r2) = (1/2,1/2). We can observe that the gap 
between the threshold and the capacity curves is between 1 and 
0.4 dB for discrepancy values A g [0,20]. This proves that 
split-extending good codes for point-to-point communications 
results in good distributed codes for cooperative communica- 
tions. 

VI. Conclusions 

We proposed a new code-design method for LDPC coded 
cooperation, which is based on a split-and-extend approach. 
First, we showed that the proposed design can be advanta- 
geously applied to existing codes, which allows for addressing 
cooperation issues for evolving standards. Subsequently, we 
introduced the concepts of threshold and capacity functions, 
and we derived density evolution equations for split-extended 
codes. Some ensemble thresholds have been presented, show- 
ing that codes optimized for point-to-point communications 
can be split-extended, so that the corresponding distributed 
codes perform close to the capacity of the relay channel for 
a wide range of discrepancy values. Optimization of split- 
extended codes will be addressed in future works. Finally, 
besides advantageous applications for cooperative transmission 
systems, the proposed design can also be used for communi- 
cation systems employing HARQ schemes with incremental 
redundancy. 
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